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Charged stripes from alternating static magnetic field 

Oleg Tchernyshyovi and Frank Wilczeki 
School of Natural Sciences, Institute for Advanced Study, Princeton, New Jersey 08540 

We motivate and perform a calculation of the energy of a cold fluid of charged fermions in the 
presence of a striped magnetic background. We flnd that a non-trivial value for the doping density 
on the walls is preferred. 



The concept of "fictitious" or "statistical" gauge fields 
has proved quite fruitful in the study of the quantum Hall 
effect. Gauge fields have also figured prominently in the- 
oretical attempts to understand other highly correlated 
electron systems in two dimensions, especially the doped 
cuprate planes famous for supporting high-temperature 
superconductivity |^|^. In the mean- field treatment of 
such theories, one is often involved in expanding around 
non-zero background values for the fictitious magnetic 
field. This background field spontaneously breaks P and 
T symmetry, and is therefore naturally associated with a 
Z2 order parameter. Unfortunately, experiments to de- 
tect macroscopic P and T violation in real materials have 
produced negative results. 

On the other hand, recently a number of experiments 
have discovered evidence for inhomogeneous electronic 
structure in the doped cuprate planes @-@|. Exist- 
ing evidence is consistent with the idea that this struc- 
ture is roughly describable in terms of antiferromagnetic 
domains separated by narrow antiphase walls, with the 
dopants concentrated on the walls, at least for small dop- 
ing. In view of this observed inhomogeneity, it is inter- 
esting and natural to consider the possibility of stripe 
structure in the context of fictitious gauge fields. Indeed, 
this concept has several advantages that are evident prior 
to any detailed calculation. The empirical difficulty of 
macroscopic P and T violation is ameliorated, and the 
Z2 nature of the order provides a raison d'etre for stable 
domain walls. 

On deeper consideration it appears that there can be, 
in the context of these ideas, a compelling connection be- 
tween doping and stripe structure. Indeed, let us suppose 
that the preferred bulk phase, realized in the absence of 
doping, prefers values of the magnetic field ±_B, with one 
sign chosen uniformly throughout the plane, at which 
the electrons just fill an integer number of Landau levels. 
As one dopes away from this preferred filling, two things 
might happen. It might be that one just has to live with 
an unfilled band. Alternatively, it might be preferable 
to produce regions with both ±_B, separated by narrow 
walls. One could have then have the filled band in bulk, 
accommodating dopants on the domain walls. 

To determine the viability of this concept, we have 
performed the following calculation. 

Consider ideal fermions in 2 dimensions in uniform 
magnetic field B. Initially, exactly N first Landau lev- 
els (n = 0, 1, . . . , — 1) are filled. The energy density 



is £0 = \NfiuJcP, where p = NeB /2Trhc is the fermion 
density and coc = eB/mc is the cyclotron frequency. 

Add or remove Ap fermions per unit area. If the mag- 
netic field stays commensurate with density, keeping N 
filled Landau levels, A£ = NfiUc Ap. Otherwise, if B 
remains fixed, 



AE = 



(N + l/2)hujcAp ifAp>0, 
(N ~l/2)hujcAp if Ap < 0. 



(1) 



In our alternative scenario, magnetic field is allowed 
to alternate between B and —B, thus producing mag- 
netic domains (Fig. |^). Landau levels are distorted in 
the vicinity of domain walls. Midgap states induced by 
the domain wall accommodate newly doped fermions (or 
holes), so that doped charges form stripes along the do- 
main walls. We wish to compare the energy of this state 
of a doped system with that in a fixed uniform magnetic 
field. 




FIG. 1. Striped arrangement of magnetic fleld. Also shown 
are the semi-classical electron orbits. Notice that states on the 
wall have a preferred direction of momentum. 

The linear density of charge on a stripe v and the stripe 
spacing A are related to the average density of added 
charge according to Ap — v jX. For small Ap, stripes are 
far apart compared to the width of the wave functions, 
and they can be considered independently of one another. 
In units of magnetic length (?ic/eS)^/^, this requires A » 
(A^- 1/2)1/2. 
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For a single domain wall of length L with = vL 
added particles, the energy is E{M) = Eo + e{L')L, where 
Eq is the energy of the undoped system without the do- 
main wall and e{v) is the energy density of a domain wall 
(per unit length). The relative energy of the system per 
unit area is 



AS 



A 



Ap- 



E{Af) - Eq 
Af 



Ap. 



(2) 



Thus the preferred linear charge density vq is found by 
minimizing e{u)/\v\, the energy per doped particle (or 
hole) . Formation of stripes is advantageous if this energy 
is lower than that in a uniform field B: 



< 



{N + l/2)hu!c for electrons, 
-(TV - l/2)huJc for holes. 



(3) 



Energy spectrum of a single domain wall: For a stripe 
along the x axis, B{y) = —Bsgn{y). It is convenient to 
work in the Landau gauge Ax = B\y\, Ay — 0, so that 
translation symmetry in x is manifest. Wavefunctions 
ipiy) exp (ikx) then satisfy the Schroedinger equation 



ny) + {\y\~kfi^{y)^2E^{y), 



(4) 



in appropriate units of length and energy. 

As a function of momentum k along the stripe, there 
are 3 simple regimes. If fc > and large one has two 
well-separated parabolic wells. Then there are doublets 
near the levels of the harmonic oscillator, with an ex- 
ponentially small splitting between the symmetric and 
antisymmetric states. If fc < and large one has a single 
wedge-like well. Energy levels start (at least) at |fc| and 
thus do not contribute to low Landau levels. Finally if 
fc = 0, one has the simple harmonic oscillator. 

Thus as k varies from -|-oo to 0, the n-th Landau level 
splits into symmetric and antisymmetric branches. The 
energy of the symmetric branch £'+(fc) is lower than n + 
1/2 for large k. At fc = 0, E+ = 2n + 1/2 and E_ = 
2n + 3/2. 

Solutions of the Schroedinger equation (^) regular at 
y — oo are parabolic cylinder functions fel 



V'(2/) 



DE-i/2{{y~k)V2) ify>0, 
±^i\y\) ify<0. 



(5) 



The spectra E±(k) are determined by a boundary con- 
dition at y = 0: ^/''(-t-O) = for symmetric states and 
'(/'(-fO) = for antisymmetric ones. All the necessary 
calculations are now fast computations for MAPLE. The 
results are shown in Fig. p[ 




FIG. 2. Energies of the symmetric and antisymmetric 
eigenstates of a stripe, as functions of k. Straight lines are 
the Fermi levels for and optimal dopant density on a stripe. 

Energy of a single stripe: Further analysis of the energy 
competition is greatly simplified by the existence of a one- 
to-one correspondence between states with and without 
a domain wall. Since fc is quantized in units ^ in both 
cases, there are no phase shifts to worry about. We will 
consider in detail the case of TV = 1 filled Landau level. 

In uniform magnetic field, all states of the lowest Lan- 
dau level in — 0, — oo < fc < +oo) are occupied. To 
obtain the same number of particles in the system with 
a single domain wall in the most economical way, we 
must partially fill the two lowest branches (n = 0), as 
follows. First fill all the symmetric states with fc > and 
all antisymmetric states with fc > 0; then slightly rear- 
range occupied and empty states within the two bands 
to achieve the state with lowest energy. The Fermi en- 
ergy and momentum are determined by the equation 
E+{-kF) = E^{kF) = Ep- Note that the Fermi level 
lies below the bottom of the symmetric n — \ band (Fig- 
urc|). 

Without doping, a stripe has a positive linear density 
of energy. At large doping levels, however, E{J\f) — Eq ^ 
—N/2 + C, as holes are doped into unperturbed Landau 
states away from the domain wall. The constant C is 
negative. This becomes evident from Figure 2 if one con- 
siders Ep just above 1/2, for then most of the remaining 
electrons are in the symmetric band, i.e., below 1/2. The 
energy per doped hole is 



E{N) - Eq 



c = C/L < 0. (6) 



Because e(i/)/|i^| > for v —0, this function has a 
minimum at a finite i^q. Moreover, condition (|^) is satis- 
fied. Thus for = 1 the stripe phase has lower energy 
than the phase with a fixed uniform magnetic field, at 
least for sufficiently low doping levels, such that stripe 
wavefunctions do not overlap. 
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FIG. 3. Energy per added hole in the stipe phase, for A'' = 1 
and 2 filled Landau levels. Both curves display shallow min- 
ima at values below the corresponding no-stripe values. 

In Figure 3, we have displayed the energy per added 
hole e(i^)/|i^| for = 1 and 2 filled Landau levels in bulk. 
For = 1, we find {vqI w 0.13 particles per magnetic 
length. 

On a lattice with flux cf) piercing an elementary plaque- 
tte, the lattice constant is a = 02. The preferred doping 
density in lattice units is therefore 



l^^ol 



0.13 particles 0.13 0^/^ particles 
magn. length lattice constant 



(7) 



If = TT, as suggested theoretically |9Hll|,|l|, we find 
|fo| ~ 0.23 particles per lattice constant. 

Discussion: In our picture, the domain walls resemble 
quantum Hall edges. They support currents, alternating 
in direction from wall to wall, in the ground state. 

Our calculation has been highly idealized, of course. 
For one thing, we have ignored any possibility of en- 
ergy intrinsically associated with the fictitious field (other 
than the implicit constraint enforcing values ±i?) . While 
this approximation is in the spirit of fictitious Chern- 
Simons fields, or of fields implementing constraints, 
which have no independent dynamics, in the absence of a 
detailed microscopic model we cannot assess its accuracy. 
Similarly, we have been shamelessly opportunistic in ma- 
neuvering back and forth between lattice and continuum. 
With all due reserve, it nevertheless seems appropriate to 
mention that the sort of model discussed here is remark- 
able and virtually unique, as far as we are aware, in pre- 
dicting a non-trivial preferred dopant density along the 
stripes. Moreover, the preferred numerical value, which 
has emerged from a dynamical calculation containing no 
disposable continuous parameters, is consistent with the 
observed one. 
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